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H ! Abstract 

We use the Paris nucleon-antinucleon optical potential for explanation of experi- 
Q\ ' mental data in the process e + e~ — ► pp near threshold. It is shown that the cross 

section and the electromagnetic form factors are very sensitive to the parameters 
■^j- . of the potential. It turns out that final-state interaction due to slightly modified 

absorptive part of the potential allows us to reproduce available experimental data. 
We also demonstrated that the cross section in nn channel is larger than that in pp 
one, and their ratio is almost energy independent up to 2.2 GeV. 
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1 Introduction 



At present, QCD can not describe quantitatively the low-energy nucleon- 
antinucleon interaction, and various phenomenological approaches have been 
suggested in order to explain numerous experimental data , see , e.g., Refs. 
[Tll2T3~ll4"ll5f6] and recent reviews [7)8] . However, parameters of the models still 
can not be extracted with a good accuracy from the experimental data [9] . 

Very recently, renewed interest in low-energy nucleon-antinucleon physics has 
been stimulated by the experimental observation of a strong enhancement of 
decay probability at low invariant mass of pp in the processes J/ty — > •ypp 
[I0|, B + -»■ K+pp and B° -> D°pp [TIT121T3] . B + -»■ n+pp and B + -> K°pp 
[13] , T — > 7pp [T3] . One of the most natural explanation of this enhancement 
is final state interaction of the proton and antiproton |16II17II18)19)20[21] . 

A similar phenomenon was observed in the investigation of the proton (an- 
tiproton) electric, Ge(Q 2 ), and magnetic, Gm(Q 2 ), form factors in the process 
e + e~ — > pp [22)23)24] , Namely, it was found that the ratio \G e(Q 2 ) / 'G m{Q 2 )\ 
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strongly depends on Q 2 = 4E 2 (in the center-of-mass frame) in the narrow 
region of the energy E near the threshold of pp production. Such strong de- 
pendence at small (E — M)/M is related to the interaction of proton and 
antiproton at large distances r >> \/M. Therefore, it is possible to take one 
of the nucleon-antinucleon interactions determined in [1.2,3,4,5 ,6J for describ- 
ing the final state interaction in the process e + e~ — > pp in order to explain the 
experimental data. In the present paper, we use the Paris nucleon-antinucleon 
optical potential V N n which has the form [4"]|5] : 

V N n = U N fj — % W N f; , (1) 

where the real part U N fj is the G-parity transform of the well established Paris 
iViV potential for the long- and medium-ranged distances (r > lfm), and some 
phenomenological part for the short distances. The absorptive part W N $ of the 
optical potential takes into account the inelastic channels of NN interaction, 
i.e. annihilation into mesons. It is essential at short distances and depends 
on the kinetic energy of the particles. Our knowledge of W N $ is essentially 
more restricted than that of U N ^. Therefore, we hope that experimental data 
for the cross section of the process e + e~ — > pp can significantly diminish the 
uncertainty in W N ft. 



2 Amplitude of the process 



Using the standard definition of the Dirac form factors Fi(Q 2 ) and F 2 {Q 2 ) [22] 
of the proton we write, in the nonrelativistic approximation, the amplitude of 
NN pair production near threshold as follows (in units e 2 /Q 2 ): 



(2E + M)F 1 {Q 2 ) + 



E 2 + 2ME 
M 



F 2 (Q 2 



EF 2 {Q 2 )-MF l {Q 2 
3M(E + M) 



3(p , -e M )p'-p% 



(2) 



where e M is a virtual photon polarization vector, and e\ is the spin-1 function of 
NN pair. Two tensor structures in Eq.(T2]) correspond to the s-wave and d-wave 
production amplitudes. The total angular momentum of the NN pair is fixed 
by a production mechanism. The functions Fi(Q 2 ) and F 2 (Q 2 ) are related in 
a standard way with the electric and the magnetic form factors of the proton 
G E = F 1 + j^pF-i and Gm — Fi + F 2 . Near threshold, it is more convenient 
to work with the Dirac form factors, because at threshold Ge = Gm- The 
amplitude (j2J) already includes effects of final state interaction. Therefore, the 
form factors in Eq.(j2J) should have a pronounced Q 2 behavior near threshold. 
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Our aim is to single out these effects. In order to do that, let us introduce the 
amplitude T\ M which has the form of Eq.pl) but with the replacement Fx — > Fx 
and F 2 — > F 2 , where the form factors F 12 do not account for the effect of 
final state interaction. So, the difference between T AAt and is i n their Q 2 
dependence due to the form factors. In T A/i , the only scale is the nucleon 
mass M (which is still too small to use perturbative QCD). Therefore, near 
threshold, where the final state interaction is important, the form factors F\^ 
are smooth functions of Q 2 and can be treated as phenomenological constants. 

Now we can construct the amplitude of NN pair production in a certain 
isospin channel / = 0, 1 using the wave function *J,V f (p') of the NN pair m 
momentum space: 



T 1 



d 3 p' 
(2^)3 p a v 

EF 2 - MF{ 
3M(E + M) 



(2E + M)F( + 



bI , E 2 + 2ME ~, 



M 



3(p' • e>' - p ,2 e 



(3) 



The Fourier transform of the function <fr pA ^(p') ; which is the wave function 

*p ( A _)t (r) of the NN pair in coordinate space, is the solution of the Schrodinger 
equation 



M 



NN i 



(4) 



where V NN is given by Eq.([T]). Note that ^(r) is the left eigenfunction of 

the bi-orthogonal set of eigenfunctions of the non-Hermitian operator H. Its 
asymptotic behavior at large distances is 



/(-)t/ 

pA I 



_ P l P r 
tpr _|_ fl_ 

r 



(5) 



For pp production, we have Tj^ = + T^, while for nn = - T^. 



2.1 Coupled channels basis 



In the presence of tensor forces, the states with angular momentum L and 
L + 2 are coupled, while the total angular momentum J is conserved. As a 
result, the final-state wave function ^(r) can be represented in the form 

= E Dt m ** aL,M \r) + E Fx JM ^Ar)Yj M (n) , (6) 

JMa JM 
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where n = r/r , 

^UM {r) = ( r)Y J-i( n ) + wf(r)Yj+\n) (a = 1 , 2) , (7) 

and 

«5(r)Y^(n) (8) 

are three independent solutions of the Schrodinger equation having total an- 
gular momentum J. In Eqs.© and jSJ) , Yj M (9,(f>) is the vector spherical 
harmonic, which is an eigenfunction of the operators L 2 , J 2 , and J z , where 
J = L + S and S = 1. In order to find the coefficients D® IJM and F[ JM in Eq. 

we substitute the asymptotics of the functions Wj 7 (V) , u>j 7 (r) and Vj(r) 
at large r in Eq. ([6]) and use Eq.([5]). Then we obtain 

Dl"M = A ^i L (JM\Lml\)Yt m {i>)A I L \ 

Lm 

F{ JM = A^ J {JM\Jml\)Y: jm {p)A^ (9) 

m 

where the coefficients A L a and Aj are related to incoming flux in the asymp- 
totics of the radial functions Uj a (r),Wj a (r) and Vj(r). Making the Fourier 
transform of the expansion Eq. ([!)]), substituting the result in Eq.([n]), and per- 
forming the integration over the angles of the vector p', we obtain 



T i 



p' 2 dp' 



1/ n »,s^t E 2 + 2ME~ r 

~(2E + M)F + Fl 

V 1 1 3M 2 



(10) 



where fj a (p') and gj a {p') are the corresponding radial parts of the wave func- 
tions in momentum space. Integrals over p' are related to the values of radial 
wave functions at r = 0. 



< a (V) = ^- 2 jv 2 fl a {v)dp , w I 1 a "(0) = - T ^- 2 Jp*g I 1 a (p)dp, (11) 

where w{ a " (0) is the second derivative at r = 0. The amplitude Eq. ljTO"]) be- 
comes 
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Using Eq.([9]) for the coefficients D, we can present the amplitude in the form 

G{<^ - v^G^|2mlA>y 2 * m (p), (13) 



1 \n 



where 



Gq 



Gl 



1/ „ , rt nr E 2 + 2ME ~ T 

-(2E + M)F + R 

3 V ; 1 3M 2 



5V2 EFj-MF{ , Ia „ 
+ 4 M(£ + M) Wl (U) 

-(2E + M)F 1+ w F 2 



5 EFj-MFj t Ia „ 
+ 2V2 M(E + M) ^ A * Wl (0) - 



(14) 



(15) 



The observable electric and magnetic form factors are expressed in terms of 
Gq and G\ in the following way 



Ge — Gq + ^^i-^, G\ { — Gq + ~^g[ . 



(16) 



3 pp production 



The proton form factors are the sum of isoscalar and isovector form factors. 
The differential cross section for pp production is 



da 



a 2 pC 



\G p M (Q 2 ) 



| 2 (1 + cos 2 i 



+ 



4M 2 



|G^(Q 2 )| 2 sin 2 £ 



(17) 



where (3 = v/c and C is the Coulomb distortion factor. We omitted here 
effects of Coulomb-nuclear interference. Using the form factors F? as free 
parameters, we fit the observed energy behavior of the cross section up to 
100 MeV of proton kinetic energy. We found that using the original Paris 
NN potential [26] it is impossible to reproduce the energy behavior of the 
cross section. The cross section falls down too steeply becoming more than 
one order of magnitude smaller than the observed one at 10 MeV of proton 
cm. kinetic energy. We verified that possible smooth dependence of the form 
factors Ff on Q in a narrow energy region near threshold (for instance, Ff oc 
l+b(Q 2 — 4M 2 )/8M 2 with b ~ 1), which is not related to final state interaction, 
does not change this result. The result is also insensitive to variation of most 
parameters of the potential. Variation of the energy independent parameters 
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is absorbed by the fitting parameters F*. However, if we modify the only 
parameter, the energy dependence of absorptive part of the triplet potential 
W N n, decreasing it by a factor of 8 -j- 10, we obtain a good fit of the cross 
section (see Fig.l). Modification of the energy dependence of the real part of 
the potential is not so important. Note that the values of the parameters, which 
are responsible for the energy dependence of the potential, are not well known. 
In the two versions of the Paris NN potential, [5] and [26], these parameters 
differ from each other up to factor 2. In order to clarify the importance of 
the energy dependence of the absorptive potential, we calculated S-matrix 
elements and corresponding phase shifts for J = 1, S = 1 at different energies. 
It turns out that the phase shifts obtained with our strong modification of 
the absorptive potential differ only slightly as compared to those obtained in 
[5]. The only noticeable modification appeared in the energy dependence of 
the parameter n directly related to absorption. In contrast, (0)| 2 , see 

Eq. ([TJ) , is much more sensitive to this modification. Therefore, the uncertainty 
in determination of the parameters of energy dependence in the absorptive part 
of the triplet potential from the scattering data is apparently larger than factor 
two, and it is necessary to use also another data. In the process e + e~ — > pp 
we have a unique situation where the quantum numbers of pp pair are fixed, 
J = S = 1. However, the absorptive part of NN effective potential is not 
universal and may be different in different scattering channels, as it takes place 
in Nijmengen potential |27j. Thus, the modification of the energy dependence 
in the absorptive part of the potential for J = S = 1 channel does not allow 
us to make any conclusions on other channels and a new fit of absorptive 
potential in these channels should be performed. 

Using the parameters obtained from the fit of the cross section, we simultane- 
ously reproduce the energy behavior of the form factors ratio \Ge/Gm\ (see 
Fig. 2). Emphasize that the form factor Ge differs from the form factor Gm 
due to the contribution of D-wave only, see Eq. ffTB"]) . Therefore, the strong 
energy dependence of the ratio \Ge/Gm\ clearly indicates the importance of 
D-wave even in the vicinity of threshold. Having obtained the amplitudes for 
the two isospins, we have calculated the cross section and the ratio \Ge/Gm\ 
for the process e + e~ — > nn. The corresponding results are shown in Figs 1,2 
by the dashed lines. It is seen that the final state interaction leads to strong 
enhancement of both quantities in nn channel as well. It is interesting that 
the cross section in nn is larger than that in pp channel, and their ratio is 
almost energy independent up to 2.2 GeV. 

Very recently, the final state interaction in e + e~ — > pp has been discussed 
in Ref. [28]. The authors of the paper have presented many arguments in 
favor of importance of final state interaction in pp production near threshold. 
Using Jiilich model for NN interaction they have calculated the contribution 
of 3 Si partial wave to the e + e~ — > pp cross section near threshold. However, 
since they have neglected the 3 D t partial wave contribution, they have not 
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been able to reproduce the ratio \Ge/Gm\ which is equal to unity in their 
approximation. 
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Fig. 1. Fit of the cross section of pp 
production (solid line). Data are from 
Ref.pl]. Dashed line: the cross section 
for nn production 



Fig. 2. Ratio of \G P E /G P M \ (solid line). 
Data are from Ref. |24j . Dashed line: 
the ratio for a neutron. 



In summary, we calculated the effects of final state interaction in the reactions 
e + e~ — > pp and e + e _ — > nn near threshold. We found that the last version 
of the Paris NN potential [26] does not reproduce the energy dependence 
of the observed cross section. A smooth dependence of the form factors F/ 
on Q in a narrow energy region near threshold, which is not related to final 
state interaction, does not change this result. Variation of the parameters 
responsible for the energy dependence of the imaginary part of the potential 
in the channel with J = S = 1, and L = 0, 2 allows us to reproduce both the 
energy dependence of the cross section an the ratio \Ge/Gm \ iorpp production. 
We obtained that the cross section in nn channel is larger than that in pp one, 
and their ratio is almost energy independent up to 2.2 GeV. 
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